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ABSTRACT 

Star-forming galaxies constitute the majority of galaxies with stellar masses > 10 10 /T 2 M Q at z ~ 2. It is thus 
critical to understand their origins, evolution, and connection to the underlying dark matter distribution. To 
this end, we identify the dark matter halos (including subhalos) that are likely to contain star-forming galaxies 
at z ~ 2 (z2SFGs) within a large dissipationless cosmological simulation and then use halo merger histories 
to follow the evolution of z2SFG descendants to z ~ 1 and z ~ 0. The evolved halos at these epochs are then 
confronted with an array of observational data in order to uncover the likely descendants of z2SFGs. Though 
the evolved halos have clustering strengths comparable to red galaxies at z ~ 1 and z ~ 0, we find that the bulk 
of z2SFGs do not evolve into red galaxies, at either epoch. This conclusion is based primarily on the fact that 
the space density of z2SFGs is much higher than that of lower redshift red galaxies, even when accounting for 
the merging of z2SFG descendants, which decreases the number density of z2SFG descendants by at most a 
factor of two by z ~ 0. Of the ~ 50% of z2SFGs that survive to z ~ 0, ~ 70% reside at the center of z ~ 
dark matter halos with M > 10 12 h~ l M Q . Halo occupation modeling of z ~ galaxies suggests that such halos 
are occupied by galaxies with M r < -20.5, implying that these z2SFGs evolve into "typical" ~ L* galaxies 
today, including our own Galaxy. The remaining ~ 30% become satellite galaxies by z ~ 0, and comparison 
to halo occupation modeling suggests that they are rather faint, with M r < — 19.5. These conclusions are at 
least a partial departure from previous work due primarily to the increased accuracy of observational data at 
Z < 1, and to higher resolution A^-body simulations that explicitly follow the evolution of dark matter subhalos, 
whose observational counterparts are likely satellite galaxies. These conclusions are qualitatively generic in 
the sense that any halo mass-selected sample of galaxies at one epoch will evolve into a more complex and 
heterogeneous sample of galaxies at a later epoch. This heterogeneity is driven largely by the fact that some 
galaxies will continue to accrete matter and form stars throughout their evolution, while others will become 
satellites and thus have their growth suppressed relative to galaxies in the field. 
Subject headings: galaxies: evolution — galaxies: halos — galaxies: high-redshift 



1. INTRODUCTION 

The past decade has witnessed an enormous increase in 
our knowledge of the high-redshift (z > 2) Universe. The 
unique rest-frame UV spectral signatures of star-forming 
galaxies have allowed observers to select galaxies at z > 
2 based solel y on optical photometry with very high ef- 
ficiency (e.g. [Steidel et all I2003L [2004). Extensive multi- 
wavelength follow-up of these objects, including optical and 
near-IR spectroscopy, and imaging at wavelengths from X-ray 
through radio, has yi elded estimates of rest-frame UV lumi- 
nosities and color s (ISteidel et alJll999t [Adelberger & Steidel 
2000trR eddv et d] |2007h , stellar masses dShaplev et alj|2005b 
Erb et al|~ 2006cl). star-formation rates ( Erbetal. 2006b: 
Reddv et al.ll2006l) . chemical abundances dPettini et alJl200lt 
Erb et al]l2006ah. and clustering strengths dAdelberg er et all 
2003Ll2005HGiayalisco & Dickinsonl200lllOuchi et aljl2005h 
Lee et al.ll2006t lHamana et alJl2006t lOuadri et al.ll2007h of a 



star-formation rate density ( Redd y~et al J 12007b at an epoch 
when the Universe was only a few billion years old. Com- 
plementary techniques tuned to the rest-frame optical prop- 
erties of galaxies at similar epochs have identified objects 
with typically redder colors and higher mass-to-light (M/L) 
ratios, but the corresponding spectroscopic confi rmation has 
only been obtained for the brightest, rare ob jects dFranx et all 
l2003llDaddi etani200l Ivan Dokkumll2006l) . 

At the same time, models connecting galaxies to the 
underlying dark matter distribution have been put forth which 
either focus on quantifying the statisti cal relation between 
galax i es and dark matter halos (e.g . iBerlind & Weinberg 



2002; Bullock et al 
Scoccimarro et all ~ 



120021 lYan et al l 12001 [Seiial l20Qi 



2001h IZhengl l2004t iKravtsov et all 



2004; IConroy et all 120061) . or pred icting the connection 



large, well-defined population of galaxies that dominate the 
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based on physical principl e s (e.g. [W hite & Frenk 1991; 
Some rville & Primackl [l999l ICole et all l2000t lHatton 1 
2003; Springe l~l.ll2001aHCroton et alfeooa Wwer 1 
2006). Essential to these latter efforts has been the vast 
increase in size and resolution of iV-body simulations, which 
has only recently allowed for the construction of detailed 
merger trees that follow not only the evolution of distinct ha- 
los, i.e. those halos not contained within any larger halo, but 
also subhalos — the halos contained within the virial radii of 
distinct halos — who se observational coun terparts are likely 
satellite galaxies (e.g. Springel et al. 2001a). Throughout we 
refer to both distinct halos and subhalos generically as halos. 
Despite impressive advances in our understanding both of the 
properties of halos themselves and of the relation between 
galaxies and halos, fundamental questions remain. 
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The origin and fate of z2SFGs is one such ques- 
tion. Herein we define z2SFGs as those galaxies se- 
lected to lie at z ~ 2 based on optical U„G1Z photom- 
etry that are brighter than 1Z = 25.5, which corresponds 
to a rest-frame UV magnitude limit at this epoch. His- 
torically, UV-selected high-redshift galaxies have been de- 
scribed as either low-mass, merger-induced starbursts or 
more massive galaxies "quiescently" forming stars. With 
early data, it wa s not possible to distinguish between these 
opposing ideas ( Lowenthal e t al.l 119971: IColes et al.l |1998[: 
Mo et all 119991; iGiavalisco & Dickinson! |200U iKolatt et all 
1999tlSomerville et al.l2001l:IWechsler et al.l2001l) . However, 
recent modeling in conjunction with more recent data on the 
abund ance and spatial clustering seems to favor t he latter pic- 
ture dConrov et al.ll2006t lAdelberger et al.ll2005l) . The distri- 
bution of stellar populations of these obj ects lends support to 
the massive, qui escent scenario as well (Shaple v" et al.ll200il 
lErb et al.ll2006cT) . 

Once identified at high redshift, it is also important to de- 
termine what type(s) of galaxies z2SFGs evolve into at later 
epochs. Most previous analyses have focused on the so- 
called Lyman-Break Galaxies, which are UV-selected star- 
forming galaxies at z ~ 3. We argue in later sections that 
these galaxies are qualitatively similar to z2SFGs and thus 
here we make no distinction between the two. Some analy- 
ses focusing on the strong observed clustering of z2SFGs fa- 
vored the scenario where z2SFGs evolved into the observed 
massive red galaxies at z ~ 1, and b y z ~ were at the 
centers of r ich groups and clusters dM o & Fukugita 119961: 
Baug h et alJll998b iGovernato et alJfl998l 120011: iBlaizot et al.1 
2004; Adelberger et alj|2005l) . though other models incorpo- 
rating the high observed comoving number density of z2SFGs 
suggested a more nuanced history ( Mo ustakas & Somervilld 
120021) . Subsequent to these efforts, there has been a vast in- 
crease in our understanding of the number densities and clus- 
tering of gal axies as a function of luminosity and color both 
at z ~ 1 fe.g.lCoil et al.ll2006tlPollo et al .1120061: iMeneux et all 
120061 ICoil et alj 12007b and z ~ (e.g. Eehavi et al.l 12005b . 
Better constraints on the statistics of galaxy populations at 
lower redshift are crucial for a robust identification of the de- 
scendants of high-redshift galaxies. These recent data at z < 1, 
in conjunction with new high-resolution simulations and dark- 
matter-halo merger trees, motivate a re-evaluation of the na- 
ture and subsequent evolution of z2SFGs. 

In the present work we focus on the fate(s) of observed 
z2SFGs. This redshift range is our focus here rather than 
Z > 3, as there has been much recent attention to this lower 
redshift window with several complementary galaxy selection 
techniques. Plus, a more extensive set of multi-wavelength 
imaging and spectroscopy exists for z ~ 2 galaxies, allowing 
the more robust estimate of a number of physical properties 
such as stellar, gas, and dyn amical masses, and extinction- 
corrected star-formation rates dErb et alj 2006c; Shaple vet alj 
l2005tlForster Schreiber et alJ2006l:lReddy et alJ2006h . which 
are important to consider along with the inferred host dark- 
matter halo properties. Such a detailed comparison of halo 
and galaxy masses and formation histories is currently not 
possible at higher redshifts. 

Under the assumption of a tight correlation between galaxy 
light and halo mass, the observed clustering of z2SFGs puts 
a constraint on the minimum halo mass hosting such galax- 
ies (fJU). Halo merger trees extracted from a cosmological 
Af-body simulation are then used to follow the halos hosting 
z2SFGs to later epochs. Comparing the clustering, satellite 



fraction, and number density of these evolved halos to ob- 
served galaxy populations at both z ~ 1 and z ~ then pro- 
vides constraints on the fates of z2SFGs (|j3]l. 

Throughout we assume a ACDM cosmology with 
(fl m , n A ,a & ) = (0.25,0.7 5,0.9), consistent with the WMAP1 
data (ISpergel et al. 2003). Where applicable, we leave results 
in terms of h, the Hubble parameter in units of 100 km s" 1 
Mpc" 1 . In several sections we discuss how our results are af- 
fected by adopti ng the cosmologic al parameters favored by 
the WMAP3 data dSpergel et alj2007l) . where the primary dif- 
ference is a lower a&. Halo masses are measured as the mass 
interior to a region with me an enclosed den sity equal to 200 
times the critical density. A lChabrierl (120031) IMF is assumed 
when quoting stellar masses. All magnitudes are quoted in 
the AB system; we omit the factor of 51og(/i) when quoting 
absolute magnitudes for brevity. 

2. PRELIMINARIES 

2.1. Connecting Galaxies to Halos 

The clustering strength of a given sample of galaxies can 
be used to estimate the minimum dark matter halo mass, 
M m i n , hosting s uch galaxies, for a specifie d ACDM cosmol- 
ogy (see e.g. IWechsler et all 1998. 2001; AdelbergeretaD 
120031: lOuadri et al.ll2007l:lGawiser et alll2007l) . This approach 
assumes that a sample of galaxies above a given luminosity 
threshold corresponds, at least approximately, to a sample of 
halos above a given mass threshold. Under this assumption, 
the threshold M m ; n is varied until the clustering of halos with 
M > M m i n matches that of the observational sample in ques- 
tion. 

While early attempts only considered distinct halos, re- 
cent work has demonstrated that this approach can be re- 
fined by including dark matter subhalos — halos that orbit 
within the potent ial wells of distinct halos — as possible sites 
for galaxies (e.g. | Colrn et alj[l 999; Kravtsov & Klypi nll 19991; 
Sprin gel et al.l 12001 at [Kravtsov et ak| [2004; Tasitsio mi et al.l 
l2004t IVale & Osu-ikel l2004l 120061: IConrov et al.ll2006h . The 
approach is as simple as before — is varied until a 
match with the observed correlation function is obtained — 
except now galaxies can reside within subhalos, and hence 
there can be multiple galaxies per distinct halo. This refine- 
ment is desirable because the observed correlation function 
requires there to b e multiple galaxies per distinct halo (e.g. 
IZehavi et al.lf2005h . 

In this model, the space density of halos with M > M m \ n 
need not match the space density of observed galaxies. This 
approach thus allows the possibility of there being fewer than 
one galaxy per halo or subhalo, which, if favored by the data, 
would indicate that z2SFGs are "on" only a fraction of the 
time dMartini & Wein berg 2001). Such a scenario may arise if 
the star-formation in these galaxies is episodic due to, for ex- 
ample, major mergers. The converse possibility, where there 
is more than one galaxy per halo, would signal an inconsis- 
tency in the model, since the model assumes that halos (both 
distinct and subhalos) can host at most one galaxy. 

Previous work has demonstrated that better agreement with 
observed small-scale clustering data can be achieved when us- 
ing the subhalo mass at the time when it is accreted on to a 
larger halo (at the "epoch of acc retion" ), rather than its p resent 
day mass dNagai & Kravtsov! 120051; IVale & OstrikeJ 120061: 
IConrov et al] l2006t IWang et al] l2006t iBerrier et all l2006h . 
The accretion epoch mass is expected to correlate more 
strongly with galaxy properties because, while a subhalo can 
experience mass-loss due to tidal stripping, the galaxy, which 
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is embedded at the very center of the halo, is much less af- 
fected by tidal evolution. For subhalos, we use this accretion- 
epoch mass herein when counting halos above M m \ n , but note 
that the dynamical times at z ~ 2 are short (~ 0.1 Gyr) and 
thus subhalos merge rapidly. Therefore, those that are identi- 
fied at z ~ 2 were only accreted rather recently. In fact, more 
than half of the subhalos identified in the simulation at z ~ 2 
above our best-fit M m ; n (see below) were accreted within the 
last simulation output, which at these epochs is ~ 0.2 Gyr. 
Because of their recent accretion, they have thus lost little 
mass due to the effects of tidal stripping; the average halo 
has lost 30% of its mass since accretion, as inferred from 
a high-resolution A^-body simulation (see below). The dif- 
ference between using accretion-epoch versus current masses 
should thus be small; we include it only for completeness. 

In sum, we vary M m ; n until the clustering of the halos with 
M > Mmin matches the clustering of the observed z2SFGs. 
We include both distinct halos and subhalos so that there can 
be multiple galaxies per distinct halo (as there are multiple 
galaxies per group/cluster). For subhalos we use the subhalo 
mass at the epoch of accretion as this should be a better proxy 
for the stellar, and hence luminous, content of the subhalo. 

As applied to z2SFGs, the approach outlined above as- 
sumes that the rest-frame UV luminosity of galaxies is tightly 
and monotonically correlated with dark matter halo mass. 
Since this does n ot appear to be the case at z ~ (e.g. 
iHeinis et al.l l2007h . such an assumption may at first glance 
appear unjustified. In fact, however, several lines of evi- 
dence suggest that this assumption is valid at z ~ 2. The most 
straightforward evidence comes from the fact that the large 
scale clustering strength of observed z2SF Gs is an increas- 
ing fu nction of rest-frame UV luminosity ( Adelberg er et al.l 
2005) at the ~ 2a level. Larger samples a t somewhat higher 
redshifts (3 < z < 5) confirm this trend (lOuchiet all 12005b 
iLee et al] |2006). Such facts are most readily understood 
if there is a correlation between UV luminosity and halo 
mass since higher mass halos are more strongly clustered 
than lower mass halos. At lower redshi ft, this general line 
of reasoning is r outinely employed (e.g. IZehavi et al.l 120051: 
IConrov et al.l2006h . and has been verified by more direct mea- 
sure ments of halo mass, such as weak gravitational lensing 
(e.g. lMandelbaum et al.ll2006l) . 

Further evidence comes from the observation of a clear cor- 
relation between the star-formation rate and baryonic mass 
(stellar mass and cold gas mass within ~ 6 kpc) for z2SFGs 
(Erb et al. 2006b, c). If the vast majority of galaxies at z ~ 2 
are forming stars vigorous ly enough to be detected as z2SFGs, 
as appears to be the case dFranx et alJ2003l:lDaddi et al.l2004t 
Ivan Dokkuml 1200 6). then a selection on rest-frame UV lu- 
minosity, which traces star-formation - modulo the effects 
of dust extinction - should correspond to a selection on 
baryonic mass. Since a strong redshift-independent corre- 
lation between baryonic and halo mass is a firm predict ion 
of any theory of galaxy formation (e.g. ICrain et al.ll2007l) . a 
UV luminosity-limited sample is approximately a halo mass- 
selected sample at these epochs. 

There is thus ample evidence that rest-frame UV luminosity 
is tightly correlated with halo mass at z ~ 2. However, the lack 
of any strong dependence of clustering on UV luminosity at 
z ~ suggests that this is not the case in the local Universe 
(e.g. IHeinis etlDl2007h . We now briefly provide a plausible 
explanation for this difference between high and low redshift. 

At high redshift, star-formation and baryonic mass are 
monotonically related primarily because most of the gas 



in halos at high redshift is cold (~ 10 4 K), a nd so it di- 
rectly contributes to the star-formation rate (e.g. iKeres et all 
2005). Only galaxies within the most massive and hence 
rarest halo s have had star-formation largely truncated at 
this epoch dFranx et alJ[2003t iDaddi et al.ll2004t Ivan Dokkuml 
2006; Oua dri et alJ 12007). At later times, infalling gas in 
massive halos is sh ock-heated to the vi r ial temperature of 
the halo (~ 10 7 K ; Iwhite & Reesl fl978t IKeres et alJ 12001 
iDekel & Birnboiml 120061: ICattaneo etaiTl2007l) . High tem- 
perature gas is much more susceptible to further heating pro- 
cesses, and thus we can assume that once gas is shock-heated, 
it remains hot forever (observations indicate that this must 
be the case for t he hot gas that perm eates groups and clus- 
ters of galaxies; iPeterson et al.l [2003b . The transition from 
predominantly low- to high-temperature gas quenches star- 
formation in high mass halos, thereby breaking the mono- 
tonic relation between star formation and baryonic (and thus 
halo) mass. These trends are generically reproduced in hydro- 
dynamic cosmolpgical simulations (e.g . iBlanton et al J 120001 
IKeres et alJl2005t ICattaneo et all 120071) . and provide a justi- 
fication of our assumption that UV light and halo mass are 
tightly correlated at high redshift, despite the absence of such 
a trend at lower redshift. 

2.2. Data atz^l 

The observed z2SFGs are selected by U„G1Z color cuts and 
required to be brighter than 1Z = 25.5. These are objects sat - 
isfying the "BX/MD" criteria of lSteidel et all (12004 120031) . 
but their comoving number density has been corrected for 
the incompleteness arising from star-forming galaxies scatter- 
ing out of the U n GTZ color selection region due to photomet- 
ric er rors (see e.g. lAdelberger et al.ll2004 120051; iReddv et al] 
120071) . This sample is thus meant to encompass all star- 
forming galaxies at z ~ 2 with 1Z < 25.5. As estimated 
from the rest-frame U V luminosity function, the number den- 
sity of these gala xies is 11 x 10~ 3 /i 3 Mpc~ 3 , with an uncer- 
tainty of ~ 10% (IReddv et alJl2007l) . This is ~ 80% higher 
than the preliminary n umber density estimates reported in 
lAdelberger etail d2005l) . 

The photometric limit in the observed 7?.-band corresponds 
to the rest-frame UV at these epochs. At z ~ 2 M*(1700A) = 
-20.2 and TZ = 25.5 corresponds to M*(1700A) = -18.6; the 
z2SFG sample thu s extends ~ 1.6 magnitudes below M* 
jReddv et alj|2007l) . This rest-frame UV limit approximately 
corresponds to an unobscured star-formation rate threshold of 
15M Q yr" 1 , given the typical z2SFG rest-frame UV exti nction 
factors of - 4 - 5 (ISteidel et al.ll2004t IReddv et alJl200l . 

The angular clustering o f these z2SFGs has been measured 
by lAdelberger et al] (120051) . who find that the observed clus- 
tering can be fit with a real-space (de-projected) correlation 
function that is a power-law, £ = (r/ ro)~ 7 , with 7 = 1 .6 ± 0. 1 
and ro = 4.2 ± 0.5 h~ l Mpc 9 . Note that while the number den- 
sity has been completene ss corrected, the c l usterin g measure- 
ments have not. That is, lAdelberger et al] ([2005) only mea- 
sured the correlation function for objects satisfying the U„G1Z 
color criteria. We are thus assuming that the objects not 
included in the clustering analysis have the same clustering 

9 The de-projection of the observed angular correlation function that 
yie lds estimates of r p and 7 is cosmology-dependent. The results reported 
by lAdel berger et al. (2005) assumed a cosmology where (f!,„, f^A, crs) = 
(0.30,0.70,0.9) — slightly different from the one adopted herein. We ex- 
pect that the change in ro induced by updating the cosmology is within la of 
the value reported above. 
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FIG. 1. — The correlation function for halos with M > 10" 4 /i _i Mq at 
Z ~ 2 (solid line) that have clustering pro perties similar to observed z2SFGs 
(hatched region; Adelberger et al. 2005). Also plotted are the correlation 
functions of those halos evolved to z. ~ 1 and z ~ 0. The correlation func- 
tion increases toward lower redshift both because of the increased clustering 
of halos and because of the higher satellite fraction at later times, which in- 
creases the weight given to higher mass (more clustered) halos. Halos with 
masses M > 10 1L2 /r 1 M Q and M > 1O 1L6 /T 1 M at z ~ 2 are included for 
comparison (lower and upper dotted lines, respectively). 

properties as those that were included. Since 1Z < 25.5 star- 
forming galaxies were missed from the color selection due 
primarily to photometric scatter ( i.e. by random pro cesses), 
this assumption is well motivated dReddv et al.ll2007l) . 

2.3. N-Body Simulations 

The properties and histories of dark matter halos are taken 
from the Millennium simulation, a large cosmological Af -body 
simulation with suffi cient resolution to i dentify both distinct 
halos and subhalos (Springel et al. 2005). Here and through- 
out distinct halos are those halos that are not contained within 
any larger virialized system while subhalos are halos whose 
centers are contained within the virial radii of larger systems. 
This classification is analogous to the distinction between cen- 
tral and satellite galaxies. The simulation box has length 
500/T 1 Mpc, particle mass m p = 8.6 x 10 8 /i _1 M Q , and was 
run with the following cosmological parameters: a% = 0.9, 
f2 m = 1 - fl\ = 0.25, h = 0.73, and n s = 1 where n s is the 
slope of the initial linear power spectrum. Halo merger trees, 
which connect halos acros s time, also exist for this simula- 
tion 10 (Springel et al. 2005), and are utilized below. In partic- 
ular, the merger trees are used to define the accretion-epoch 
mass for subhalos as the mass when the subhalo was last con- 
sidered a distinct halo. Working with a simulation of this 
size can be computationally expensive; for this reason we thus 
make use of a 150 3 /T 3 Mpc 3 region of the simulation. This 
size is sufficient for our purposes; for example, making use of 
a 250 3 /T 3 Mpc 3 region at z ~ 2 has a negligible effect on our 
conclusions at that epoch. 

In order to investigate the sensitivity of our results to cos- 
mological parameters, we performed two additional Af-body 
simulations. One used the same cosmological parameters 
as the Millennium run, and the other used cosmological pa- 
rameters from the WMAP3 results (tt m = 0.239, Q A = 0.761, 

10 The merger trees and halo catalogs are publically available and can be 
found here: http : / /www . mpa-garching . mpg . de /millennium. 



Q h = 0.04166, h = 0.73, n s = 0.953, erg = 0.756: ISpergefeTall 
120071) . Each simulation contained 512 3 particles in a box 
of length 141.3 h~ l Mpc, resulting in a particle mass of 
1.46 x 10 9 /i _1 M Q for the Millennium cosmology and 1.39 x 
I0 9 h- 1 M Q fo r the WMAP3 cosmology. The GADGET-2 
N-body code (ISpringel et al.ll2001bl 120051) with a Plummer- 
equivalent softening length of 7.9h~ l kpc was used for both 
simulations. 

In these additional simulations, we identified distinct halos 
with a friends-of-friends algorithm 11 using a linking length of 
0.2 in units of the mean interparticle separation. The halos of 
interest contained at least 70 particles. No subhalo detection 
was performed on these simulations. We define halo mass 
as the mass within a sphere centered on the minimum poten- 
tial energy halo particle (after removing unbound particles) 
and enclosing a density equal to 200 times the critical density. 
This definition closely matches the definition of M cr it,200 use d 
in the Millennium simulation (which first removed subhalo 
particles before identifying the halo particle with minimum 
potential energy) and elsewhere in this paper. 

The descendant of a z ~ 2 halo was identified by finding 
the halo at a later time that contained the plurality of the par- 
ticles comprising the z ~ 2 halo. However, the position of the 
descendant was taken to be the later-time center-of-mass of 
the particles comprising the z ~ 2 halo. We explicitly verified 
that our simulation with the Millennium cosmology produced 
consistent distinct halo mass functions, clustering as a func- 
tion of halo mass, and halo descendant clustering as a function 
of halo mass as the Millennium simulation. 

3. RESULTS 
3.1. The z2SFG '-Halo Connection 

The clustering of halos extracted from the simulation at 
z ~ 2 is shown in Figure Q] for our best-fit M m ; n threshold, 
along with the observed z2SFG clustering. As discussed in 
H2.3\ here and throughout "halos" refers to both distinct ha- 
los and subhalos. The clustering of halos with M m i„ ±0.2 
dex from the best-fit M m [ n is included for comparison. It is 
clear from the figure thatM m i n = 10 11 A h~ l M Q provides a good 
match to the observed clustering of z2SFGs. The slope of the 
halo correlation function on scales 1 < r < I0h~ l Mpc is -1 .5 
which is consistent at the lcr level with the slope inferred from 
observations (-1 .6). This figure also contains the clustering of 
the z ~ 2 halos evolved to z ~ 1 and z ~ 0, based on the Mil- 
lennium simulation merger trees. The space density of z ~ 2 
halos with M > M min is 7.5 x 10" 3 /? 3 Mpc" 3 . 

The uncertainties of the z2SFG clustering data translates 
into an uncertainty on M m ; n of ~ ±0.2 dex (as determined 
by eye; see Figure [TJ. This uncertainty affects both the in- 
ferred number densities and the clustering strength of the cor- 
responding host halos. ForM m i n = 10 1L2 h~ l M Q the number 
density is 13 x 10 _3 /i 3 Mpc" 3 while for M min = 10 11 ' 6 h~ l M Q 
it is 4.3 x 10~ 3 /i 3 Mpc -3 . This range brackets the observed 
number density of z2SFGs (1 1 x 10~ 3 /i 3 Mpc" 3 ) and indicates 
that, within the uncertainty, every halo above M m ; n contains 
one z2SFG. The uncertainty in M m \ n is explicitly incorporated 
into our uncertainty in the clustering of z2SFG halos and their 
descendants in the following sections. 

In reality the dark matter halo occupation function of 
z2SFGs need not be a step function at M m \ n (zero galaxies 

" The algorithm we use has been made freely avail- 
able by the University of Washington HPCC group: 

http : / /www-hpec .astro . Washington . edu/tools/f of .html. 
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per halo/subhalo below M mm and one above it). For exam- 
ple, scatter between galaxy UV luminosity and halo mass at 
z ~ 2 will result in a more gradual rise of the occupation func- 
tion from zero to one around M mm . In order to understand 
the qualitative effect of scatter on our inferred number den- 
sity of halos that match the clustering of z2SFGs, we have run 
a series of halo occupation models that were c onstrained to 
match the observed correlation function (see e.g. Tin ker et al.l 
2006, for details). In these models both M mm and the amount 
of scatter between mass and light were left as free parame- 
ters. The resulting number density of halos with clustering 
matching the observed z2SFGs varies by ~ 20% compared to 
the number density of halos when setting the scatter to zero 
(which is our default model herein). The effect of scatter is 
thus insignificant for our purposes. 

Using ro to constrain the minimum dark matter halo 
mass of high-redshift ga laxies is not a new technique 
(e.g. I Wechsler et all 1 19981) . Previously, Adelberg er et al. l 
(120051) used the GIF-ACDM simulation (Kauffma nn et al.l 
1999) to constrain the minimum mass, finding that M m i n ~ 
10 1L8 h~ l Mq provided a good fit to the z ~ 2 data (which is 
the same data used herein). The difference between this value 
and ours (10 1L4 /z -1 M Q ) is due to the updated ft,,, and more 
physically motivated transfer function in the simulations we 
use. 12 At first glance this may be worrisome, since the cos- 
mology used herein is alre ady somewhat out- dated insofar as 
the new WMAP3 results (Spergel et al.ll2007l) favor a lower 
normalization of the power spectrum (erg = 0.76) than what 
we assume. To explore the dependence of our conclusions on 
cosmological para meters, we have made use of simulations, 
described in £12.31 with the updated WMAP3 parameters and 
find that the best-fit M mm decreases by 0.3 - 0.4 dex. This 
lower M m i n results in a 35 - 70% increase in the abundance 
of z2SFG host halos — well within the uncertainties asso- 
ciated with matching halos to galaxies using clustering, as 
discussed above. The differences between the WMAP1 and 
WMAP3 cosmological parameters thus have a small impact 
on the halo-z2SFG connection. 

In sum, the observed clustering strength and comoving 
number density of z2SFGs combine to suggest that every dark 
matter halo and subhalo with mass > lO rL4 /! _1 M is host to 
roughly one z2SFG. The subhalo fraction among these halos 
implies that 1 1% of z2SFGs are satellites. 

Techniques complementary to the optical color-selection 
method used to identify z2SFGs have identified a population 
of massive, red galaxies at z ~ 2 with low space density (i.e. 
objects with n ~ 10~ 4 /i 3 M pc" 3 ; iFranx et alJ2003t iDaddi et al l 
l2004t IvanD okkum 2006). A significant fraction of these ob- 
jects may be missed by the z2SFG UG7Z criteria — espe- 
cially if they have little or no ongoing star-formation. How- 
ever, since the observed z2SFG sample has a space density of 
n ~ 10~ 2 /z 3 Mpc" 3 , our results are not sensitive to this possibly 
distinct population of massive red galaxies because such ob- 
jects are comparatively rare and thus have little impact on the 
number density and clustering of the overall set of z2SFGs. 

3.1.1. Quiescent versus Collisional z2SFGs 

In the previous section we found that the number density 
of halos above M m [ n was consistent with the observed num- 
ber density of z2SFGs, given the uncertainty in M m j n . Thus, 

12 The GIF simulations used an analytic fitting function for the power 
spectrum transfer function (Bond & Efstathiou 1984) while the Millennium 
simul ation utilized the more accurate CMBFAST code (Seliak & Zaldarriaga 
119961) to generate it. 



within the uncertainties, all halos above M m { n contain one 
z2SFG. Such a scenario arises naturally if z2SFGs are qui- 
escently forming stars, since in this picture every sufficiently 
massive halo should have the same conditions necessary for 
star-formation (e.g. a sufficient supply of cold gas). Here and 
throughout, "quiescent" refers to an approximately continu- 
ous star-formation history, as opposed to an episodic one. 

However, it is less clear that a collisional starburst scenario 
(where star formation is merger-induced and thus episodic) 
would be consistent with this result, since the majority of 
halos, at any mass scale, are not undergoing violent merg- 
ers frequently enough. For example, based on merger trees 
in the Millennium simulation, only ~ 20% of halos above 
Mrmn = 10 nA h~ l Mq at z ~ 2 have had a merger with mass 
ratio < 3 : 1 in the past gigayear. If elevated star formation 
occurs for such mergers but not for mergers at larger mass 
ratio s, as suggested by controlled hydrodynamic simulations 
(e.g. iMihos & Hernquist|[l994t ICox etafl 120071) . then there 
are simply too few such mergers in halos with M > M m [ n to ex- 
plain the abundance of observed z2SFGs under th e collisional 
starburst scenario. Previou s modeling efforts (Kol att et al.l 
119991; IWechsler et"aT] 1200 ll) concluded that the collisional 
starburst scenario was viable because they included collisions 
with very high mass ratios (10 : 1 and greater). It now ap- 
pears, however, that such merger ratios will not result in the 
enhanced star-formation ( Cox et af]l2007h necessary to make 
them detected as z2SFGs. Moreover, the small scatter in the 
observed star-formation rate-stellar mass relation at z ~ 1 sug- 
gests that by this epoch starbursts a re not the dominan t mode 
of star formation in most galaxies dNoeske et al. 2007). 

Of course, to accomodate the low rate of roughly equal 
mass mergers, one could lower M m j n until the number den- 
sity of halos with low mass-ratio mergers (< 3 : 1) matched 
that of the z2SFGs. However, lowering M m \ n in order to in- 
corporate more colliding halos would also lower the large 
scale clustering stength, and would quickly produce disagree- 
ment with the observed clustering of z2SFGs (note that any 
"merger bias", in the sense that colliding halos are more clus- 
tered than other halos of the sa me mass, is probably small; see 
Scannapieco & Thacker 2003). Models that rely on collisions 
as the impetus for star-formation would thus have a very dif- 
ficult time simultaneously matching the observed clustering 
and space density of observed z2SFGs. 

Some degree of caution is in order when attempting to 
discriminate between these two scenarios for the origin of 
z2SFGs, as the accretion and merging rate of halos is signifi- 
cantly higher at high redshift compared to z ~ 0, and thus the 
dichotomy between quiescent and collisional star-formation 
may not be applicable at high redshift. As noted above, major 
mergers, which are the types of mergers thought to fuel colli- 
sional starbursts, are still sufficiently rare to rule that channel 
out as a major contributor to the star-formation rate. But mi- 
nor mergers of ratio ^10:1 are far more common at higher 
redshift, and far fewer halos are truly quiescent, at least in 
terms of their accretion of dark matter. For example, the ha- 
los above M m ; n at z ~ 2 in the Millennium simulation have 
a mass doubling time of ~ 1 Gyr; only the rarest and most 
massive halos have s uch vigorous accretion properties today 
(IWechsler et aUl2002h . 

3.1.2. The Baryon Budget 

Here we briefly compare the available baryon reservoir in- 
ferred from our results to estimates of the gas+stellar mass in 
a subsample of observed z2SFGs. 
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The average halo mass for the z2SFGs is 10 11 • 8 A~ 1 Af ; this 
can be converted into an average baryonic mass by assuming 
that each halo contai ns the universal baryon fraction, fb = 0.17 
( Sperg eTet al.l 12007). Under this assumption, the average 
baryonic mass for the z2SFG sample is 1.1 x 10" h~ l Mq. For 
a fraction of the z2SFGs modeled in the present work there 
exist stellar and gas mass estimates (lErb et al.ll2006d) . The 
observed gas ma sses were estimated by utilizing the empirical 
iKennicutTd 1 9981) law that relates star-formation rate densities 
to cold gas densities. The star-formation rate was estimated 
from the Ha emission line and its surface density was esti- 
mated from the ~ 6 kpc extent of the emission. Stellar masses 
were estimated from stellar population synthesis modeling 
with a lChabriedd2003l) IMF. The average observed gas+stellar 
masses of these objects is 5.8 x 10 10 M Q , for h = 0.7. 

In order to compare our expected average baryon fraction to 
this gas+stellar mass estimate we convert our mass to h = 0.7 
units, and find that the observed gas+stellar mass estimate is 
roughly a factor of three smaller than the average baryonic 
mass inferred from the halos. Based on the assumption that 
the Kennicutt law holds at z ~ 2, and thus that the estimated 
gas masses are accurate, our results indicate that on average a 
third of the baryons in the halo are in the central ~ 6 kpc. Fur- 
thermore, while > 50% of the observed baryons in z2SF Gs 
appear to have been converted into stars (lErb et al.ll2006cllal) . 
the inferred baryonic mass from the average z2SFG halo mass 
indicates that only ~ 20% of the total available baryons have 
been converted into stars. Moreover, halos are constantly be- 
ing fed new gas as they accrete gas-filled dark matter halos. 
There is thus a sufficient supply of gas to sustain star forma- 
tion for many gigayears, if the gas is able to cool. As we 
will see below, this feature is important because many of the 
z2SFG descendants are galaxies that are still forming stars at 
Z ~ 1 and z ~ 0. 

3.2. z2SFG Evolution t o z ~ 1 and z ~ 

The identification of the halos likely containing z2SFGs 
is the first step toward understanding the evolution of these 
galaxies to z ~ 1 and z ~ 0, where there are well-defined 
samples su ch as the Deep Ex tragalactic Evolutionary Probe 
2 (DEEP2; iDavis et alJl2003h. th e VIMOS VLT Deep Sur- 
vey (VVDS: iLe Fevre et al.ll2005l). the 2dF Galaxy Redshift 
Survey (2dFGRS: IColless et alj 120011). and the Sloan Digi- 
tal Sky Survey (SDSS:lAdelman-McCa rthv et alj|2006l) . The 
next step is to evolve these halos forward in time using the 
halo merger trees from the Millennium simulation. Note that 
the merger trees follow the evolution of subhalos until they 
have been destroyed, i.e., until the halo finder no longer iden- 
tifies the subhalo as a bound clump of particles. 

The primary ambiguity in this next step is how we treat 
z2SFG subhalos that, according to the merger trees, have 
merged/disrupted at some later time. These cases are ambigu- 
ous because the Millennium simulation may, like any sim- 
ulation, artificially destroy a subhalo , either because of r es- 
olution limits (e.g. iMoore et alJll999t iKlypin et alJll999h or 
because the simulation does not include the effects of baryon 
condensation, which may make a subhalo more resilient to 
disruption in the real universe (although the latter effect is ex- 
pected to be small for the regimes of interest here; see e.g. 
INagai & Kravtsovll2005l: IWeinberg et al.ll2006l) . We incorpo- 
rate this uncertainty by computing results in this section for 
two cases that should bracket the range of possibilities. In 
the first case we assume that the simulation is correct and 
that when a subhalo merges, so does the galaxy embedded 
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FIG. 2. — Relationship between clustering strength (ro) and sample number 
density (n) for observed galaxies and the descendants of the halos hosting 
z2SFGs at z ~ 1 (top panel) and z ~ (bottom panel). The data at z ~ 1 
are for samples defined above various magnitude thresholds (from —19.5 to 
-20.5 in half magnitude steps for the overall sample, and from —19 to —21 
for the color-defined samples; Coil et al. 2006, 2007), while at z ~ they 
are defined for ma gnitude bins (in one magnitude intervals from —19 to —22 
Zehavi et al. 2005). We have also included data from the "all" sample at z ~ 
for magnitude threshold samples (diamonds) in ord er to show the differences 
between binned and threshold samples (see i|3.2. ll for details). For the halos, 
eiTor bars encompass t he un certainty due both to M mm and the merger vs. no 
merger scenarios (see £|3.2l for details). The location of the symbols for the 
halos indicates the values for the merger scenario withM mm = \Q UA h- l MQ 
(large open boxes). z2SFG descendant halos with a satellite fraction lowered 
to match observations are also included (large crossed boxes). 



within it. In the second case we assume that none of the 
subhalos actually merges. In this case, if a subhalo merges 
within a distinct halo according to the merger tree, we place 
the satellite within the di stinct halo with a po sition specified 
by an NFW distribution (Nav arro et al-l fl 997) appropriate for 
the background dark matter 13 . These will be referred to as 
the merger and no-mer ger scenarios below. For reference, the 
conclusions reached in Adelb erger et alj d2005l) . namely that 
z2SFGs evolve into massive red galaxies by z ~ 0, were based 
on the no-merger scenario. 

In what follows we focus on three constraints that will 
help discriminate between possible evolutionary histories 
of z2SFGs; these are the space density of galaxies, their 
large scale (1 < r < 10/T 1 Mpc) clustering strength, and 
the fraction of galaxies that are satellites. These three 
constraints, when combined, strongly disfavor any scenario 
where z2SFGs evolve into a single class of objects at lower 

13 In fact, the galaxies associated with disrupte d subhalos will l ikely be 
more centrally concentrated than the dark matter I Sales et al. 2007). How- 
ever, the radial distribution of galaxies within the halo has a negligible ef- 
fect on the large-scale clustering strength and leaves the satellite fraction and 
number density of the galaxies unchanged. This uncertainty thus does not 
impact our analysis. 
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redshift (i.e. red, blue, central satellite). As we discuss the 
evolution of z2SFGs to lower redshifts, it is worth remember- 
ing that z2SFGs are defined according to a luminosity-limit 
in the rest-frame UV, while samples at lower redshifts are de- 
fined according to increasingly redder rest-frame bands (B- 
band at z <~ 1 and r-band at z ~ 0). At increasingly shorter 
wavelengths, the light emitted by galaxies is increasingly 
dominated by young stars and hence recent episodes of star- 
formation, while longer wavelengths are dominated by older 
stars and hence probe the total stellar mass of a galaxy. At first 
glance then, connecting galaxies selected by star-formation 
at z ~ 2 to galaxies selected more closely by stellar mass at 
z ~ would seem to be a daunting task. However, with the 
assignment of z2SFGs to dark matter halos, connecting these 
galaxies to their lower redshift counterparts becomes simpler 
thanks to the halo merger trees, which provide a clear connec- 
tion between high and low redshift. 

The space densities, large-scale clustering strengths, and 
satellite fractions of z2SFG descendants and various observed 
samples are plotted in Figures |2] and [3] The uncertainties on 
ro and / sat for the z2SFG descendant halos reflect the un- 
certainty in Mmin and the uncertain fate of satellites within 
merged subhalos. In contrast, the uncertainty in the number 
densities of z2SFG descendant halos reflects the error in the 
observed z2SFG number density, with the additional uncer- 
tainty due to mergers. In other words, the observed z2SFG 
number density is multiplied by the fraction of z2SFG halos 
that survive to z ~ 1 and z ~ in order to deduce the number 
density of z2SFG descendants at these epochs. The uncer- 
tainty in the halo number density due to the uncertainty in 
M m i n thus is not included in the uncertainty in the observed 
z2SFG descendants. After discussing these constraints, we 
compare the distribution of z2SFG descendant halo masses 
with the halo masses of observed galaxies at lower redshifts, 
as inferred from halo occupation modeling. The most basic 
observational constraint, the abundance of galaxies of various 
types at multiple epochs, is considered first. 

3.2.1. Number Density 

In the no-merger scenario, the comoving number density 
of z2SFGs is constant with time, by construction. Including 
the possibility of z2SFG mergers (as determined by the merg- 
ing/disruption of their associated subhalos), results in only 
mildly more evolution than the no-merger case: by z ~ 1 the 
number density of z2SFG descendant halos is 76% of its z ~ 2 
value and by z ~ it is 53% of that at z ~ 2. Recall that 
our merger prescription associates a z2SFG merger with the 
merging/disruption of its associated subhalo, which is likely 
an over-estimate of the true merging because of simulation 
resolution effects. Nonetheless, if we use the fraction of sur- 
viving halos as a proxy for the galaxy number density at lower 
redshift, more than half of the z2SFGs exist as independent 
entities at z ~ 0. 

Figures [2] and [3] plot along the x-axis the number density 
of z2SFG descendants at z ~ 1 and z ~ 0. Figure [2] also 
plots along the x-axis the abundances of observed galaxies 
as a function of luminosity and color for absolute Z?-band 
magnitude thresholds at z ~ 1 (in half magnitude steps from 
-19.5 to -20.5 for the overall sample, and from -19 to -21 
for the color-defined samples) and for absolute r-band mag- 
nitude bins at z ~ (in one magnitude intervals from -19 
to -22) based on the luminosit y func tions at these epoch s 
(IWillmer et alJ 120061 iBell et alJl2004t iBlanton et al.l l2003h . 
For reference, at z ~ 1 M* = -20.6 and at z ~ M* = -20.4. 
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FIG. 3. — Relationship between satellite fraction and sample number den- 
sity for the z2SFG halo descendants (large boxes) and various observations 
(circles; symbols are as in Figure ff), at both z ~ 1 (top panel) and z ~ 
(bottom panel). Also included are z2SFG descendant halos with a satellite 
fraction measured from the Millennium simulation (large open boxes) and 
lowered to match observations (large crossed boxes). For the halos, error bars 
encompass the unce rtainty due both to M ra ; n and the merger vs. no merger 
scenarios (see ^13.21 for details). The location of the symbols for the halos 
indicates the values for the merger scenario with M m \ n = 10 l" 1 Mq. 



At both epochs, blue and red galaxies are separated accord- 
ing to the observed bimodal distr ibution of optical colors 
dBaldrv et alj|2004t IBell et al.ll2004l) . We have also included 
data at z ~ for overall magnitude threshold samples in order 
to demonstrate the small difference between using magnitude 
bins and thresholds. 14 

The mild evolution of the z2SFG descendant halo num- 
ber density already places strong constraints on the possible 
progeny of z2SFGs. It has been proposed that red galaxies 
are the descendants of z2SFGs. However, there are simply 
far too many z2SFG descendants compared to the observed 
number density of z ~ 1 red galaxies for this to be the case, 
even when including rather faint (Mb = -19.5) red galaxies. 
Since the faint-end slope of the re d galaxy luminosity func- 
tion at z ~ 1 is shallow (a = -0.5: IWillmer et al.ll2006h . the 
number density of all red galaxies to a limiting magnitude 
of M B = -18.0 (4.4 x 10~ 3 /i 3 Mpc" 3 ) still does not equal the 
abundance of z2SFG descendant halos. The abundance of an 
overall (i.e. no cut on color) sample of galaxies brighter than 
Mb = -20.0 at z ~ 1 is much more similar to the abundance of 
z2SFG descendant halos. At z ~ the conclusions are sim- 
ilar. The observed sample at z ~ with the number density 
closest to the z2SFG descendant number density is the overall 
-21 < M r < -20 sample. 

These conclusions, based solely on the abundances of halos 

14 Magnitude bins, rather than thresholds, are used at z ~ because ro is 
only available in magnitude bins for the color-defined samples. 
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and galaxies, suggest that at most a fraction of z2SFGs have 
evolved onto the red sequence by either z ~ 1 or z ~ (see also 
Moustakas & Somerville 2002; Gilli et alj2 007. who reached 
similar conclusions). The z2SFGs also do not appear to 
evolve exclusively into the most luminous (> L*) population 
of galaxies at either epoch, as the number density of such 
galaxies is much lower than the descendant halos of z2SFGs. 
Rather, the abundances of descendant z2SFG halos is most 
similar to typical ~ L* galaxies at both z ~ 1 and z ~ 0. 

Previous work by Adelberge r et al.l d2005l) came to the con- 
clusion that by z ~ 1 z2SFG descendants had largely evolved 
into red galaxies. This conclusion was based in part on an 
earlier estim ate of the observe d number density of red galax- 
ies at z ~ 1 dChen et al.l l2003). which was a fact or of several 
larger than current, more accurate estimates dWillmer et al.l 
2006). The number density of red galaxies at z ~ was sim- 
ilarly somewhat hig her than the more recent data used herein 
(Zeha vTet al.ll2005 | ). The lower number density of observed 
z2SFGs used in [Adelberger et al.1 (l2005h . compared to what 
we adopt from more recent data, also made it appear in that 
work that z2SFGs evolved into red galaxies by z ~ 0. Our 
use of more accurate and up-to-date number densities at all 
epochs is a significant reason w hy our results are a dep arture 
from the conclusions reached by Adelberger et al.l d2005l) . 

3.2.2. Clustering Strength 

Figure [2] plots the large-scale clustering strength, ro, as a 
function of number density for both observed galaxy sam- 
ples and z2SFG descendant halos at z ~ 1 {top panel) and 
z ~ {bottom panel). The data samples include both over- 
all and color-defined samples, for a range of luminosities de- 
scribe d above. The cl uster ing data fo r these samples come 
from |Coileiai1 (I2Q06J) a nd ICoil et all d2007l) for z ~ 1 and 
from lZehavi et al.1 d2005) for z ~ 0. For the halos, the power- 
law fit to the correlation function is computed by solving si- 
multaneously for ro and 7 over the range 1 < r < I0h~ l Mpc. 
The uncertainty on ro for the halos is dominated by the uncer- 
tainty on M m i„ at z ~ 2. 

From the figure it is clear that the clustering of z2SFG de- 
scendant halos is comparable to that of the most luminous 
red galaxies, at both epochs. This point has been used as ev- 
idence that z2SFGs ev olve into red galaxies by z ~ 1 (e.g. 
Adelberge r et al.1 120051) . However, as discussed in the pre- 
vious section, the constraints from the abundance of z2SFG 
descendants strongly disfavors this scenario. Moreover, ro is 
sensitive to the number of satellite galaxies, and thus should 
not be used as a constraint without also considering the ob- 
served constraints on the satellite fraction. 

3.2.3. Satellite Fraction 

We now turn to a discussion of the satellite fractions of 
z2SFG halos and their descendants, and of observed galaxies 
as inferred from halo occupation modeling. For z2SFG host 
dark matter halos, the satellite (i.e. subhalo) fraction at z ~ 2 
is 11%, while at z~ 1 andz-Oit is 24-39% and 31 -59%, 
respectively. The lower and upper bounds in satellite fraction 
depend primarily on the assumption of merging or no merging 
of the satellites within destroyed subhalos. Figure[3]compares 
these satellite fractions t o satellite fractio ns derived for vari- 
ous data sets at z ~ 1 (IZheng et al.ll2007l) 15 and z ~ 0. The 

15 Zhena et al. 1 2007) define a halo to be the mass enclosed within a region 
200 times the mean density of the Universe (rather than the critical density as 
is used herein), and so their satellite fractions will be somewhat larger than 
those calculated using our halo definition. 
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FIG. 4. — The average number of z2SFG descendants per distinct halo at 
z ~ as a function of distinct halo mass for the overall descendant sam- 
ple (filled circles) and for the satellites only (diamonds). These results 
are compared to the average number of galaxies per halo as inferred from 
halo occupation modeling, for all galaxies with magnitudes M r < —20.5 and 
M r < — 19.5, and for satellites only. 



derived satellite fractions were estimated by simultaneously 
fitting the observed number densities and projected two-point 
correlation functions of the galaxy samples with halo occu pa- 
tion models (see e.g. Tinker et al. 2005; Zheng et al. 20Q3, for 
details). The satellite fractions for the color-defined samples 
at z ~ are reported here for the first time. 

From this plot it is clear that the z2SFG descendant halos 
have higher satellite fractions than any overall luminosity- 
limited sample at either z ~ 1 or z ~ (except perhaps the 
faintest sample at z ~ 0). Only the fainter (-21 <M r < -19) 
red galaxy samples at z ~ have comparable satellite frac- 
tions. On the other hand, these faint red galaxies have large- 
scale clustering strengths that are lower than the z2SFG de- 
scendant halos (Figure [3J and they are thus not the likely de- 
scendants of z2SFGs. The satellite fraction of z2SFG descen- 
dants is higher than that of a typical mass-selected sample of 
halos because, while the sample is mass-selected at z ~ 2, it is 
not mass-selected at lower redshifts due to the different evo- 
lutionary histories of subhalos and distinct halos (described in 
more detail in £|3.2.51 l. 

In order to better illustrate how the satellite fraction and 
clustering strengths are related, we lower the satellite frac- 
tions of the z2SFG descendant halos at z ~ 1 and z ~ so 
that they agree with the observed satellite fractions of ~ L* 
galaxies at both epochs. We do this by removing a fraction 
of the oldest subhalos (i.e. those subhalos that accreted at the 
highest redshifts) such that the resulting halo satellite frac- 
tions are 17% at z ~ 1 and 15% at z ~ 0. In other words, at 
z~0we have removed 60% of available z2SFG descendant 
satellites. As seen in Figure|2j it is clear that the clustering of 
the z2SFG descendant halos with lower satellite fractions is in 
much better agreement with the clustering of observed over- 
all ~ L* galaxies at both epochs. This result is not strongly 
affected by removing the oldest subhalos as opposed to a ran- 
dom set; we remove the oldest because it is these subhalos 
that are most likely to have faded significantly in luminosity, 
and thus may plausibly not be included in a sample of ~ L* 
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FIG. 5. — Mass function for z2SFG distinct halo descendants {solid line), 
z2SFG subhalo descendants (dotted line), and all distinct halos (dashed line), 
at z ~ 1 (top panel) and z ~ (bottom panel), not including subhalos. For 
subhalos, the mass refers to the mass it had at the epoch of accretion. Note 
that for masses > 10 12 ft - ' Mq the halos that contained z2SFGs constitute the 
vast majority of all halos above this mass at z ~ 1 and z ~ 0. 



galaxies at lower redshift. Moreover, reducing the satellite 
fraction to zero results in ro values of 4.2 and 4.9/r 1 Mpc at 
Z ~ 1 and z ~ 0, respectively. These numbers bracket the ob- 
served values near ~ L*. The implications of this exercise are 
discussed more fully in the sections that follow. 

3.2.4. Insights from Halo Occupation Modeling 

Halo occupation modeling has emerged as a power- 
ful tool to analyze o bserva tio nal data (e.g. ISeliakl 
Berli nd & Weinberg] l2002t IScoccimarro et al 
Bullock et all l2002t iKravtsov et aTT 



2000; 



2001 

2004; Zehav i et al 



2005t ItMCTetal] l2005t ICooravl 120061; ITinker et al.l 120071: 
Zhen g et al.l 120071: Ivan d en Bosc h et al.l 20071) . The model 
quantifies the statistical distribution of galaxies within dark 
matter halos, specifying the probability that a halo of mass 
M contains N galaxies of a particular type. A given halo 
occupation distribution maps nearly uniquely onto a two- 
point correlation function, so the occupation function can be 
constrained by clustering measurements. In the context of 
this model, "halos" refers only to distinct halos, and not the 
subhalos within them. 

The function (N{M)) specifies the mean number of galax- 
ies in a halo of mass M and is typically split into two terms, 
one describing the number of central galaxies, (N c ), and the 
other describing satellite galaxies, (N s ) (the dependence on 
halo mass is implicit in these functions). When constraining 
(N(M)), the approach taken is to assume that (N c ) rises from 
zero to one at some mass scale, with a rapidity constrained 
by the data. The number of central galaxies never rises above 
one by definition. The number of satellite galaxies is taken to 
be a power-law in halo mass, (N s ) oc M a , where a is usually 



found to be near unity. This form is well -motivated both by 
results from hvdrod vnamic (|Zheng et all 20051) and dissipa- 
tionless simulations ( Kr avtsov et al. EOOIT and observations 
(e.g. iLin et al.ll2004l) . Details of the hal o occupation result s 
presented in this se ction can be found in ITinker et al.l (120051) . 
ITinker et al.l (120061) . and lTinker etaTI (120071) . although the fits 
presented herein are new because they utilize a halo mass def- 
inition that is identical to the mass definition used throughout 
the rest of this paper - the mass enclosed within a region that 
is 200 times the critical density of the Universe. 

Figure 2] plots the (N{M)) that best matches the abundance 
and clustering of observed galaxies from the SDSS survey 
with magnitudes M r < -20.5 and M r < -19.5 at z ~ 0. The 
figure shows explicitly the contribution of satellites to the full 
(N(M)). The contribution due to centrals can be inferred by 
subtracting the satellite contribution from the total. These re- 
sults are compared to the average number of z2SFG descen- 
dants per halo, where the contribution due to z2SFG descen- 
dants that are satellites is also included. 

It is clear that (N(M)) for the z2SFG descendants that are 
satellites by z ~ agrees very well with that of observed 
satellite galaxies with M r < -19.5. Moreover, the fact that 
the overall (N{M)) for z2SFG descendants drops rapidly to 
zero at roughly the same mass scale as observed galaxies with 
M r < -20.5 indicates that the z2SFGs that are central galax- 
ies by z ~ correspond closely to galaxies with M r < -20.5. 
This last fact follows because the region where (N{M)) < 1 is 
dominated by central galaxies. Scatter in the relation between 
halo mass and UV luminosity at z ~ 2, as discussed in &!3.1I 
will produce a softer roll off in (N(M)) at lower masses for the 
z2SFG descendants than what is shown in Figure|4] The effect 
of scatter is thus to associate a fraction of z2SFG descendants 
that are centrals with ga laxie s fainter than M r = -20.5. How- 
evever, as discussed in Q2.1\ the scatter is not expected to be 
substantial because of the observed UV-luminosity-dependent 
clustering observed at z ~ 2. 

This comparison vividly demonstrates that z2SFG descen- 
dants that are satellites by z ~ evolve into observed faint, 
M r < -19.5, galaxies, while those descendants that remain 
distinct halos (i.e. central galaxies) evolve into more lumi- 
nous, M r < -20.5, galaxies. Echoing the conclusions drawn 
from previous sections, it is clear from this figure that z2SFG 
descendants do not generically evolve into a single class of 
galaxies at later epochs. Rather, the evolution of z2SFG de- 
scendants that become satellites is qualitatively different from 
those that become central galaxies. 

The ratio between (N s ) for galaxies brighter than M r = 
-20.5 andM f = -19.5 is ~ 30%, indicating that approximately 
30% of galaxies brighter than M r = -19.5 are also brighter 
than M r = -20.5. In §323\ we found that keeping only 40% 
of the total number of z2SFG descendants that are satellites 
yielded a sample of z2SFG descendants that had both clus- 
tering strengths, satellite fractions, and abundances in good 
agreement with the observed sample of galaxies at z ~ with 
M r < -20. The fact that these fractions are similar suggests 
a physical motivation for reducing the satellite fraction of the 
z2SFG descendants in order to produce better agreement with 
an observed z ~ sample defined with respect to a luminos- 
ity cut. A large fraction of z2SFG descendants that become 
satellites are faint, and so the overall sample of z2SFG de- 
scendants, which includes these faint satellites, will not com- 
pare to any observed sample of galaxies defined according to 
a simple luminosity cut. 
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3.2.5. Descendant Halo Masses 

Finally, we turn to a discussion of the mass distribution of 
halos hosting z2SFG descendants at lower redshifts. Figure 
[5] plots the mass function of distinct z2SFG descendant ha- 
los at z ~ 1 and z ~ and compares to the full mass function 
of distinct halos from the simulation. It is clear that z2SFG 
descendant halos with mass > IQ 12 Ii~ 1 Mq constitute the vast 
majority of all halos in that mass range at z ~ 1 and z ~ 0. In 
other words, virtually all halos today with mass > 10 12 /r'M© 
contained at least one z2SFG. In light of this fact, it thus ap- 
pears that the descendants of z2SFGs are not only the galaxies 
at the centers of massive halos but rather constitute a wide va- 
riety of objects in a variety of environments. In fact, by far 
the most common place to find the descendants of z2SFGs at 
z ~ is in ~ 10 12 /! _I M Q halos. 

It is furthermore clear from the figure that the average 
mass of distinct z2SFG descendant halos is growing mod- 
estly with time. At z ~ 2 the average distinct halo mass is 
10 1L8 /i _1 M o , atz~ 1 it has increased to 1O 121 /T 1 M and by 
z ~ it is 10 12 ' 3 /i _1 M Q . Note that these average masses are 
not weighted by the number of satellites/subhalos within each 
distinct halo. 

Figure [5] also includes the mass function of z2SFG descen- 
dants that are subhalos by z ~ 1 and z ~ 0. Here the sub- 
halo mass is measured at the epoch of accretion. Aside from 
the obvious fact that the subhalos constitute only a fraction 
(/sat) of the total halo population, the mass function of these 
subhalos clearly peaks at a lower mass than for the distinct 
descendant halos. This is not surprising given the fact, for 
example, that the average redshift of accretion for subhalos 
identified at z ~ is z ~ 0.5. The mass that these subhalos 
had at the epoch of accretion is on average only 55% of the 
average mass of z ~ z2SFG descendant halos that are not 
subhalos. Thus, in the ~ 3.5 h Gyr since z ~ 0.5, z2SFG 
descendants residing at the centers of distinct halos contin- 
ued to grow, both in dark matter, gas, and presumably in stars 
(via star formation), while those in subhalos did not grow in 
dark matter or gas, and were likely subject to one or more 
gas-starvation processes, thereby further halting the growth 
of the galaxy embedded within the subhalo. This difference 
between z2SFG descendants that are subhalos and those that 
are distinct halos has significant implications for the galaxies 
likely embedded within them. 

3.2.6. Sensitivity to Cosmological Parameters 

As discussed in iQ.ll we have made use of an A^-body sim- 
ulation run with cosmological param eters advocated by the 
WMAP3 results ( Spergel et al.l 120071) in order to ass ess the 
sensitivity of our results to these new parameters. In i Q.ll we 
found that the halos in the WMAP3 cosmology with clustering 
most similar to observed z2SFGs have M m m lower by 0.3-0.4 
dex compared to the WMAP1 cosmology. However, once we 
have identified a population of halos with comparable clus- 
tering to the observed population in both cosmologies, their 
subsequent evolution to lower redshift is almost identical. In 
particular, we have followed the evolution of distinct halos 
(not subhalos) in the WMAP3 cosmology and find that ro of 
the descendants is 5- 10% lower compared to the WMAP1 
cosmology at both z ~ 1 and z ~ 0. Furthermore, the fraction 
of halos that merged away by z ~ 1 and z ~ differs by less 
than 5% for the two cosmologies. Note that since we use the 
observed number density of z2SFGs multiplied by the fraction 
of halos that survive to later epochs, rather than the number 



density of halos themselves, our quoted number densities of 
descendants are insensitive to cosmological parameters (see 
^.2\ for more details). The differences between these two 
cosmologies thus has a negligible impact on our conclusions 
regarding the descendants of z2SFGs. 

4. DISCUSSION 

In the preceding section we found that the observed abun- 
dance and clustering of z2SFGs implied that every halo at 
z ~ 2 above a mass of 10 11,4 /i _1 Mq contains one z2SFG. 
By z ~ 0, the overall population of halos that once hosted 
z2SFGs, and survived to z ~ 0, has a comoving number den- 
sity comparable to that of ~ L* galaxies (M r < -20; while 
M* = -20.4), a clustering strength comparable to luminous 
(Mr < — 21) red galaxies, and a satellite fraction comparable 
to ~ L* or fainter red galaxies. We found that a subsample 
of z2SFG descendant halos where 60% of the satellites were 
removed resulted in much better agreement with the cluster- 
ing strength, satellite fraction, and number density of ~ L* 
galaxies at z ~ 0. Comparison to halo occupation model- 
ing of observed galaxy samples revealed that the z2SFG de- 
scendants that are central galaxies most closely resemble ob- 
served galaxies brighter than M r ~ -20.5 while the satellite 
properties compare more favorably to galaxies brighter than 
M r ~ -19.5. The average redshift at which z2SFG descen- 
dant satellites at z ~ were accreted onto their parent halo 
was z 0.5, and the mass of the halo associated with these 
satellites when they were accreted was only 55% of the mass 
of distinct halos at z ~ 0, on average. This indicates that the 
growth of subhalos (and the galaxies within them) was signif- 
icantly less than that of distinct halos. 

These results have straightforward implications regarding 
the connection between z2SFGs and their descendants at z ~ 
0. First, it is abundantly clear that z2SFGs do not evolve into 
any single sample of galaxies defined either according to a lu- 
minosity cut or a luminosity and color cut. In particular, they 
do not evolve exclusively into massive red galaxies, either by 
z ~ 1 or z ~ 0. Rather, the subsequent evolution of z2SFGs 
is more nuanced. Merger trees extracted from simulations in- 
dicate that the number density of z2SFG descendants at z ~ 
is ~ 50% of the number density of z2SFGs because some de- 
scendants merge together between z ~ 2 and z ~ 0. This is an 
upper bound on the evolution of the number density because 
halos in our simulations can merge or disrupt without the sub- 
sequent merger or disruption of the galaxy within it due to 
simulation resolution effects or possible baryonic effects. Of 
the remaining descendants, ~ 70% are centrals and ~ 30% 
are satellites. The central galaxies have properties similar to 
observed galaxies with M r < -20.5. In particular, they reside 
predominantly in halos with mass M > 10 12 /! _1 Mq. 

An anecdotal consequ ence is that our ow n Galaxy, with 
halo mass ~ 10 12 h~ l M & (Kl ypin et alJ20 02). was once likely 
a z2SFG. This conclusion is in accord with various proper- 
ties of the bulge of our Galaxy. For example, the bulge has 
a stellar mass of ~ 1O 1O M that is thought to have formed 
over a short period of time (< 1 Gyr) r oughly 10 Gyr ago 
dZoccali et al.ll2003t [perreras et al. 2003). These facts imply 
that the star-formation rate during the formation of the bulge 
was > IOMq yr" 1 . The epoch of formation of the bulge and 
its high star-formation rate imply that it woul d have likely 
been detected as a z2SFG (see also Pettini 2006, who reached 
similar conclusions). 

In contrast, the z2SFG descendants that become satellites 
by z ~ are much more similar to a fainter sample of ob- 
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served galaxies (M r < -19.5). This is not surprising in light 
of the different histories of satellite compared to central galax- 
ies as inferred from the histories of their dark matter halos. 
Since satellites fell into their parent halo on average at z ~ 0.5 
(3.5 h~ l Gyr ago), their growth was significantly retarded rel- 
ative to central galaxies. Satellites grew less not only be- 
cause they could not accrete new gas but also because the 
gas they possessed at infall was likely prohibited from cool- 
ing to form stars (or removed altogether) due to one or more 
quenching processes, such as gas strangulation, harassment, 
or ram-pressure stripping. 

Not all satellites are required to be so faint; our results indi- 
cate that ~ 40% are brighter than M r ~ -20.5. It is plausible 
that the ~ 60% required to have -19.5 < M r < -20.5 were the 
satellites accreted at the earliest epochs, because such satel- 
lites would have the least amount of growth between z ~ 2 and 
z ~ 0. Simulated star-formatio n histories coupled with a stel- 
lar population synthesis code ( jBruzual & Charlolll2003l) con- 
firm that galaxies with properties similar to z2SFGs can easily 
evolve into galaxies with M r ~ -19.5 if their star-formation 
is truncated at z ~ 0.5 - 1 . For example, a galaxy at z ~ 2 
with stellar mass 10 I0 /i~ 2 M Q and a star-formation rate of 10 
Mq yr" 1 that steadily declines to zero at z ~ 1 fades from 
M r = -21 .5 to M r = -19.2 by z - 0. 

Similar conclusions hold for the descendant halos at z ~ 1 . 
In particular, from analysis of the Millennium simulation 
merger trees we find that ^25% of z2SFG descendants have 
merged away by this time. Of the remaining halos, ~ 25% 
are satellites and ~ 75% are central galaxies. If half of the 
satellites are fainter than ~ L* , then the remaining descendant 
halos that are satellites and centrals have a number density, 
clustering strength, and satellite fraction comparable to ob- 
served ~ L* galaxies at z ~ 1 . 

There has been some debate in the literature over whether 
or not the optical color selection technique used to select 
z2SFGs misses a significant p opulation of redd e r (i.e. older 
and/o r dustier) gala xies (" e.g. iFranx et al.ll2003t iDaddi et all 
12004 Ivan Dokkuml 120061: ISmail et alj|2002l) . Under our as- 
sumption of a tight correlation between rest-frame UV lumi- 
nosity and halo mass, we have demonstrated that at z ~ 2 
roughly every halo with mass above 10 11 ' 4 h~ l M Q contains 
one z2SFG. It would be rather surprising if each halo also 
contained another comparably massive galaxy not detectable 
with the optical color selection technique, if only by anal- 
ogy with lower redshift, where the vast major ity of ~ L* 
galax ies live alone in their dark matter halos (e.g. Zhe ng et ail 
120071) . Thus, while there may well be a population of massive 
(M* > \0 n h~ 2 Mq), red galaxies with a correspondingly low 
space density (~ 10~ 4 /z 3 Mpc~ 3 ), that is missed with the op- 
tical color selection techn ique (Franx et al. 2003; Dad di et all 
120041; Ivan Dokkuml2006l) . it is hard to imagine that this could 
be the case at lower galaxy masses. Since our conclusions are 
most sensitive to lower mass galaxies, they are not influenced 
by this possible distinct population of more massive, less nu- 
merous, red galaxies. 

There are two basic facts about dark matter halos and their 
evolution that, when combined, provide both a qualitative un- 
derstanding of our results and more general insights into the 
evolution of populations of galaxies across time. The first 
fact is that, to with a factor of ~ 2, the space density of a 
set of halos above a mass threshold does not change from 
high redshift to the present. For example, of the halos with 
M > l0 UA h- l M G at z ~ 2, 50% survive to z ~ 0. For halos 
more massive than 10 11 s h~ l M Q and 10 124 h~ l M Q at z ~ 2, 



the fraction that survive to z ~ increases to 60% and 70%, 
respectively. The second fact is that, by and large, the ranking 
of halos by mass is preserved as the Universe evolves. For ex- 
ample, the most massive halos at z ~ 2 evolve to the most mas- 
sive halos today. If, as we have argued, the z2SFGs are an at 
least approximately halo mass selected sample, then these two 
basic facts make less surprising the conclusion that z2SFGs 
evolve largely into a roughly halo mass selected sample of 
abundant galaxies at z ~ 0. As described above, it is primarily 
the uncertain evolution of satellite galaxies that complicates 
this simple description. 

These facts can also shed light on the likely descendants of 
the massive red galaxies discussed above. If such galaxies are 
in the most massive halos (as suggested by their high stellar 
mass and clustering strength; Ouadri et al. 2007), then it is 
these galaxies that will evolve, largely intact, to the centers 
of rich groups and clusters by z ~ 0. Better constraints on 
the clustering and abundances of these observed massive red 
galaxies are required to make more detailed statements. 

Finally, these basic facts can connect our discussion herein 
to the fates of star-forming galaxies at z > 2. Indeed, much 
of the previous high-redshift modeling efforts discussed in 
the Introduction focused on star-forming galaxies at z ~ 3, 
the so-called Lyman-Break Galaxies (LBGs). If we apply 
our methods to mat ch the observed cluster ing of LBGs (ro = 
4.0±0.6/T 1 Mpc; lAdelberger et all 120051) with a set of ha- 
los in the Millennium simulation above some minimum mass, 
we find a best-fit M min = 10 111±0 - 2 /T 1 M Q . Analysis of the 
z2SFG halo merger trees indicates that every z2SFG halo con- 
tains one or more of these LBG halos in their history and ap- 
proximately 70% of L BGs evolve into z2SFG s (in agreement 
with earlier work; e.g. Adelberger et al. 2005). While the link 
between LBGs and z2SFGs appears strong, there are impor- 
tant subtleties one must keep in mind when comparing LBGs 
to z2SFGs. In particular, the clustering of these two popula- 
tions has been measured for all objects above a common ap- 
parent magnitude limit (1Z = 25.5), and thus the LBGs, being 
more distant, are an intrinsically more luminous sample com- 
pared to the z2SFGs. However, the luminosity f unctions of 
these two classes of galaxies are nearly the same (Red dv et alj 
12007b . suggesting that if the two samples were defined with 
respect to the same absolute luminosity limit, their evolution- 
ary bond would be even stronger. Nonetheless, the clear link 
between LBGs and z2SFGs implies that LBGs share the same 
varied fates as those we have outlined for z2SFGs. 

The most significant uncertainty in our analysis is the clus- 
tering strength of z2SFGs and the subsequent merging of their 
descendants. A more quantitative understanding of the fate of 
these z2SFGs will thus require a more accurate accounting 
of their clustering strengths, on the observational side, and a 
more sophisticated treatment of z2SFG mergers on the theo- 
retical side. In addition, clustering measurements on smaller 
scales (< IMpc) and over a range of rest-frame luminosities 
will allow a more detailed modeling of the connection be- 
tween z2SFGs and halos at z ~ 2 than has been presented here 
and will provide a critical test of our assumption of a tight cor- 
relation between rest-frame UV luminosity and halo mass at 
z~2. 

With more accurate data over the interval < z < 2, the 
goal will ultimately be to model the galaxy-halo connection 
as a function of galaxy properties such as color and luminos- 
ity, or star-formation rate and stellar mass, in narrow redshift 
intervals. The galaxy-halo connections at higher redshift can 
then be evolved into the connections at lower redshift with 



12 



CONROY ET AL. 



the aid of halo merger trees. The merger trees are the key to 
this approach, as they put the static (fixed redshift) halo oc- 
cupation models into motion. Such an approach has only re- 
cently become feasible, due both to the vastly increased data 
at these epochs and to a converging set of cosmological pa- 
rameters that removes significant uncertainties i n the proper- 
ties and histories of dark matter halos (see e .g. IZheng et~aT] 
l2007t IWhite etalJl2007t IConrov etalJl2007h . A simplified 
version of this approach has been explored in the present 
work. A more detailed version will be pursued elsewhere. 
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